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Abstract. The diabatic approach to collective nuclear motion is reformulated in the local-density appro- 
ximation in order to treat the normal modes of a spherical nuclear droplet analytically. In a first application 
the adiabatic isoscalar modes are studied and results for the eigenvalues of compressional (bulk) and pure 
surface modes are presented as function of density and temperature inside the droplet, as well as for 
different mass numbers and for soft and stiff equations of state. We find that the region of bulk instabilities 
(spinodal regime) is substantially smaller for nuclear droplets than for infinite nuclear matter. For small 
densities below 30% of normal nuclear matter density and for temperatures below 5 MeV all relevant 
bulk modes become unstable with the same growth rates. The surface modes have a larger spinodal 
region, reaching out to densities and temperatures way beyond the spinodal line for bulk instabilities. 
Essential experimental features of multifragmentation, like fragmentation temperatures and fragment- 
mass distributions (in particular the power-law behavior) are consistent with the instability properties 
of an expanding nuclear droplet, and hence with a dynamical fragmentation process within the spinodal 
regime of bulk and surface modes (spinodal decomposition). 

PACS. 21.60.Ev, 21.65.+f, 25.70.Mn 
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1 Introduction and summary 

In bombardments of nuclei with light and heavy ions, 
as well as by absorption of pions and antiprotons in nu- 
clei, transient systems are produced at excitation energies 
around 10 MeV/u, which decay into several fragments of 
intermediate size. Increasing interest |Q,||J3[^ in this mul- 
tifragmentation reaction has been initiated by the obser- 
vation J5| that the fragment-mass distribution is propor- 
tional to A~'^ (power law with a « 2.6) indicating that 
the process may be related to the critical point of the 
liquid-gas phase transition of nuclear matter. Indeed, sim- 
ulations within molecular dynamics |6|j7|,|8|| and mean-field 
approaches |9|, [lO| , [ll|Jl2| show that initial compression in 
central heavy-ion collisions or pure excitation of a nu- 
cleus by abrasion of nucleons or by absorption of light 
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ions, pions or antiprotons cause the system to expand and 
to break up into pieces in a low-density regime. Power- 
law fragment-mass distributions have also been obtained 
from dynamical nucleation in a thermodynamically un- 
stable nucleonic system [^. However, as indicated by the 
great success of statistical models [ p^|p5| , p^ , p7t , the ob- 
servables in multifragmentation seem to be mainly deter- 
mined by the available phase space at the point of freeze 
out |18|| . This conclusion remains valid even if the time 
evolution (expansion) of the statistically emitting source 
is taken into account p9[ . 

With the discovery of a caloric curve in projectile frag- 
mentation by the ALADIN collaboration [^ , the nuclear 
liquid-gas phase transition as a mediator between the heat- 
ed expanding nuclear matter and its decay into pieces 
moved again into the center of discussions. The growth 
of density instabilities within the spinodal region inside 
the phase coexisting regime, known as spinodal decompo- 
sition |gl|], has early been considered as a possible mech- 
anism of multifragmentation p2,p3|,p4l. Extending Lan- 
dau's Fermi-liquid theory to finite temperatures Heisel- 
berg, Pethick and Ravenhall |^ have investigated the 
stability and instability of hot and dilute infinite nuclear 
matter. Subsequent studies have been concerned with de- 
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tailed prope rties of nuclear matter [E6l and finite droplets 
P7| , p8y29| , |30| in and around the spinodal regime. Quantal 
effects, as studied within constrained RPA [Q, do not 
alter the adiabatic instabilities as calculated from a fluid- 
dynamical description p7| and simulations p9| within a 
stochastic mean field approach. Some experimental evi- 
dence for spinodal decomposition has been reported by 
Rivet et al. [|l|. 

Our model for examining the stability and instability 
of a hot and dilute nuclear droplet is based on the diabatic 
approach to dissipative large-amplitude collective motion 
P^ , ^3|j3^ and, by the use of a local energy-density approx- 
imation, is similar in spirit to the Fermi liquid-drop model 
of Kolomietz and Shlomo |Q. However, instead of using 
the fluid-dynamical approach with the boundary of zero 
total pressure on the fixed surface as in |^,^ , we consider 
a general displacement field, which is defined by an expan- 
sion of the displacement potential in terms of multipoles, 
and include Coulomb interactions. This expansion with a 
suitable boundary condition for the complete basis set al- 
lows the analytical evaluation of collective mass and stiff- 
ness tensors within a consistent harmonic approximation. 
The set of eigenvalue equations couple modes with differ- 
ent number of nodes in the radial function of the displace- 
ment field. In contradistinction to ||27| , where dispersion 
relations of the energy as function of radial wave number 
for different multipoles have been obtained in the adia- 
batic limit (no deformations of the local Fermi sphere), 
we determine the orthogonal eigenmodes of the droplet as 
function of the relaxation time r for the decay of deforma- 
tions of the local Fermi sphere, i.e. continuously from the 
adiabatic to the diabatic limit. Furthermore we study also 
pure surface modes and compare the instability properties 
for soft and stiff equations of state. In this way we extend 
the exploratory studies P7|j2q , |29| , p0| of spinodal instabil- 
ities of finite droplets to a more systematic analysis of 
stability and instability by examining the eigenmodes of a 
nuclear droplet as functions of densities and temperatures, 
relevant in the multifragmentation process. 

In the following, which forms part I of a series of pub- 
lications, we present detailed results on isoscalar bulk and 
surface modes with fixed ratio of neutron and proton den- 
sities in the adiabatic limit (very fast equilibration, r — > 
0). Since r ex T^^, cf. Q, this corresponds to the high- 
temperature limit. We have studied these adiabatic iso- 
scalar modes in detail, because they are related to ther- 
modynamics and to many studies performed in the past. 
In subsequent papers we will present the modifications 
due to dissipation (0 < t < oo) and extend the study to 
the inclusion of isovector modes, which to our knowledge 
have been studied so far only for infinite nuclear matter 
ll^,!!]!!] and for a lattice gas [|o) . 

The results on adiabatic isoscalar bulk instabilities (pre- 
liminary calculations have been presented earlier [^) are 
summerized as follows. 

— As compared to infinite nuclear matter the spinodal 
region for compressional (bulk) instabilities shrinks to 
smaller densities and temperatures with Tcrit — 6 MeV 
for a soft equation of state (EOS). This effect is due to 



an increase of stability, which results from the Weizsac- 
ker term ex (Vgf)^ in the energy density and the finite 
wave lengths of density fluctuations in the droplet. The 
observed fragmentation temperatures of about 5 MeV 
(cf. pO|) are consistent with spinodal decomposition 
after expansion. Typical values for the growth rates are 
7^5 MeV corresponding to growth times h/j w 40 
fm/c. 

— Effects from Coulomb interactions on the bulk insta- 
bilities are negligible, although, as shown in E2|, they 
are important at later stages of the multifragmentation 
process, where statistical models apply. 

— As compared to the soft EOS, a stiff EOS yields a 
larger spinodal region with Tdit = 8 MeV. Typical 
values for the growth rates of instabilities are larger 
by almost a factor 2. 

— With decreasing density and temperature the modes 
with the lowest multipolarities and no radial node be- 
come unstable first. 

— At densities below 0.3go (with qq = 0.16 fm"'^) the in- 
stability growth rates for different multipolarities {I = 
2, 3, 4, 5) and number of nodes (n = 0, 1, 2, 3) are prac- 
tically equal. This property can yield a power-law be- 
havior A""' with fj Ri 2.0 of the fragment-mass distri- 
bution in agreement with experimental observations, 
cf. [O and is not related to the critical point. Ad- 
ditional components from evaporation at lower excita- 
tion energies and non-equilibrium coalescence at higher 
excitation energies are considered responsible for the 
observed larger values of a. 

For finite nuclear droplets surface modes are important in 
addition to the compressional modes. Indeed, pure surface 
modes show some interesting features. 

— The instability region of pure surface modes extends to 
larger densities up to about the spinodal line of infinite 
nuclear matter and to large temperatures. 

— In general the growth times are smaller by half an order 
of magnitude as compared to the typical values for bulk 
instabilities. 

— In the stable region surface modes are slow, such that 
deformations initiated in the excitation process will 
persist during expansion and clustering. 

— The surface instability is dominated by quadrupole de- 
formation. Fission at high excitation energies is ex- 
pected to take this path through a low density stage, 
where the fission barrier vanishes. This is a novel mech- 
anism of fission and may be even faster than the stan- 
dard path across the barrier at equilibrium density. 

In the following section we introduce the collective 
model, which serves as the basis for our studies of sta- 
bility and instability of nuclear droplets. In section y we 
describe the application to adiabatic isoscalar modes of 
compressional and pure surface deformations. Results on 
stability and instability of expanding nuclear droplets are 
presented in section 0. Relations to multifragmentation 
are discussed in section ^ Details of derivations and eval- 
uations are given in the Appendix. 
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2 The collective model 

We apply the diabatic approach to dissipative coUective 
motion ||32| , p3|j34[| . In contradistinction to adiabatic mod- 
els, dynamical distortions of the Fermi sphere are included 
and dissipation due to two-body collisions give rise to a 
term in the collective equation of motion, which is non- 
local in time (non-markovian friction). 



2.1 Reminder of the diabatic approach 

If we introduce a set g = {gx} of collective variables g\{t) 
in the irrotational velocity field 



v(r, t) ^ VVF(r, g, g) = ^ gx ^wxir, g) 



(1) 



of collective motion, we can rigorously define (stationary) 
diabatic single-particle states |<^a(g)) by 



d 
dg\ 



\Ms)) 



-AiVwx)-V + V-iVwx)}\Ms)) (2) 



i.e. by an infinitesimal unitary transformation, which de- 
forms the wave function in accordance with the velocity 
field (cf. Appendix |A|). We search for the solution \^{t)) 
of the many-body Schrodinger equation by expanding 



\Ht)) 



L^ 



c.(i)|'Z'.(g,g,t)) 



(3) 



in terms of Slater determinants \i'v)i which are built from 
the boosted diabatic single-particle states 



IV'o) 



: exp 



mwM^(g,g)- I di'ea(i') 


} \K{^)) 




(4) 



moving with the velocity field v. Here, e^ and ?7^^r denote 
the diabatic single-particle energies and the nucleon mass, 
respectively. 

From the variational principle 



■'Cu,g\ 



'dt{Wmh^^-HMt))=0 (5) 



with the many-body Hamiltonian _ff, one derives equa- 
tions of motion for the expansion coefficients c^{t) and 
collective variables gx{t). The coupled equations for the 
expansion coefficients describe the mixing of states (es- 
sentially due to residual interactions), and hence the equi- 
libration process within the intrinsic degrees of freedom. 
We approximate this equilibration by a relaxation equa- 
tion for the diabatic single-particle occupation probabili- 
ties 



dna{t) 1 



T-aiS, fJ.,T)} 



(6) 



where t and n^, respectively, denote the relaxation time 
and the equilibrium values of the occupation probabilities. 



The chemical potential /i(i) and the temperature T{t) in 
the Fermi distribution for 11^ are determined by the con- 
servation of mass and excitation energy in the relaxation 
process. The relaxation time r has been estimated for nor- 
mal nuclear matter by several groups, cf. p4| , and given 
as 

T rj 

h e* 



(7) 



with values for 77 between 0.15 and 0.30 depending on the 
magnitude of the nucleon- nucleon cross-section in medium. 
The collective equation of motion, as obtained from 
the variation (||) with respect to qx, is given by Q 

d ^-^ . 1 v^ dBx'X" 

T7 2^ i^AA' 5A' - ^ 2^ 



AA' 



A'A" 



d^A 



-9X'9\" = ^A , 



(8) 



where Bxx' denotes the collective mass tensor 

Bxx' ^mN d^rg{r){\/wx) ■ (Vu;a') 
for the velocity field (1^) and the density 
g{r) =^?i„|^„(r)|2 . 

a 

The force on the r.h.s. of (ph is determined by 

dSais) 



^^ = -E' 



dg\ 



■na{t) 



(9) 



(10) 



(11) 



in the diabatic representation. Note that the diabatic single- 
particle energies are smooth functions of g and that the 
force is defined for fixed tt-q,, i.e. for constant (single-particle) 
entropy. 

Inserting in (O) the formal integral of (0) 



= n.,m- CM''^ern"^^"/<n 



na{t) = na{t) - I At -— 



dt' 



(12) 



with the initial value na{ta) — na{ta), we obtain 

Fx{g,t) = Fx[^,t)+F'x{^,t), (13) 



F.(g,t) = -^ 



^ga(g) 

dgx 



(14) 



F'xig,t)^Yl 



dgx 



di'^e-A'd*"/^(t"). (15) 
dt' ^ ' 



The equilibrium force is determined by the equilibrium 
distribution rTc((t) for the time-dependent temperature and 
chemical potential. For temperatures (T ss 1 MeV) large 
compared to the residual single-particle coupling, i.e. dif- 
ference between the total mean field potential (adiabatic 
potential) and the diabatic potential, Fx is identical to the 
adiabatic force, because the sum of neighboring energies 
is the same in both cases (invariance of trace). The elasto- 
plastic force F'^ describes giant elastic vibrations with fre- 
quency u) for r 3> cij^^ and Markov dissipation (ordinary 
friction) for r <C oj^^ . Whereas l arge-amp litude collective 
motion has been considered in |32,p^,pJ|, we specify the 
following discussion to small amplitudes. 
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2.2 Small-amplitude modes 

For small arnplitudes around some equilibrium point go=0 
we keep in m) only terms linear in g\, and hence we have 

Y. Bxx'gy + C'^y I At'e-^'-'"^'^gx'(t') + Cxx' gy = 

A' -^io 

(16) 
with B\x' = -BAA'(go), T = T{to) and the stiffness tensors 



Caa' - Yl 



d^Sa _ dsa dria 



dgxdgx' " dgx dgy 



g=0 






(17) 



(18) 



g=o 



Inserting gx oa exp(— zwt), we obtain from (16) the eigen- 
value equation 



Gll,{Lu) = -Bxx'u;'' + C',y 



OJ 



uj + i/i 



+ Cxy = (19) 



for t — to ^ T, where Gaa' i'-^) denotes the Fourier trans- 
form of the Green function for (l6f). In general, i.e. for 
finite r-values, there are three complex solutions to this 
equation. For the two elastic (isentropic) limits, t = 
(adiabatic limit) and r — > cxd (diabatic limit) only two, 
either real (w^ > 0) or imaginary (w^ < 0) frequencies 
survive. In the diabatic limit the stiffness tensor is given 

by 



C. 



XX' 



C[ 



XX' 



Cxx' — 






dgxdgx 



(20) 



g=o 



whereas it is Caa' in the adiabatic limit. As mentioned 
below (|l5| ) the adiabatic stiffness tensor can be expressed 
for temperatures T > 1 MeV by 



where tgg and v'^p^p denote the kinetic-energy part of 
the single-particle energy and the antisymmetrized two- 
body interaction matrix element, respectively. The total 
intrinsic kinetic energy -Ekin(g) and the interaction en- 
ergy £'int(g) can be expressed by Skyrme functionals e[g\ 
for the energy density, where g denotes the deformed den- 
sity due to collective motion. Furthermore, we introduce a 
phenomenological surface energy tension es(g). Then the 
stiffness tensors are defined by 



Caa' — 



Caa' — 



d^ 



dgxdgx 



dgxdgx 



d^ref 



dfesid) 



d're[g\+ df es{g) 



g=o 



g=o 



(23) 



(24) 



in the local density approximation, where e differs from 
the adiabatic value e by the additional deformation of the 
local Fermi sphere in the diabatic limit. Note that the 
derivatives in (|2^) and (|24| ) have to be taken at constant 
entropy, i.e. for fixed occupation probabilities in accor- 
dance with (^ and (|2l|). 



2.4 Energy densities and surface tension 

We write the Skyrme energy-density functional pql as a 
sum 



<^m = fr + ey + £w + ec 



(25) 



of terms related to the intrinsic kinetic energy (cr), includ- 
ing the momentum-dependent part of the interaction, the 
nuclear interaction of homogeneous systems (ev), the con- 
tribution from inhomogeneity (Weizsacker term ew) and 
the Coulomb interaction (ec)- Explicitly these terms read 



6. = el") 



Ap) 



(26) 



Caa' ~ Y 



d'e-. 



(21) 



.95A%A'/g=0 

where a denote the adiabatic single-particle states. 

2.3 Local-density approximation 

In applying the formulation to the modes of a spherical 
nuclear droplet of homogeneous density g and temperature 
T inside the sharp surface, we perform a local density 
approximation. The mass tensor (0) is already given in an 
appropriate form. For the stiffness tensors (|2(]| ) and ( pi|) 
we make use of the relation 



Edsg 
„ dgx 



EUtgg 



dgx 



d 1 
dgx 2 



7 . '"al3a 



pn/sug 



al3 



d 
dgx 



(^kin(g)+^int(g)) 



(22) 



ev = 2^" 



d3fc fi'e -,,, 

[2-k)-^ 2m* 



(1 + 2^o)^ - [xo + ^){Qi + qI) 



12 



hQ° 



(l + :^a:3)?^-(x3 + -)(e^ + ^) 



ew = ^{-iti-t2){^gf 
Id 

+ ^{ih+t2)[{Vgn? + {Vgpf 



^ 2./ Ir-r'l 



(27) 



(28) 



(29) 



(30) 
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where g = 2 is the spin degeneracy factor, fi{k) denote the 
local momentum distributions (normalized to density) for 
neutrons (i — n) and protons {i — p) inside the droplet 
and m* the corresponding effective masses given by 

— = + [{tl+t2)Q+-{t2-h)Q^]. 31 

ml TOjv 2 

Note that the r-dependence in eqs. ([2q-P9|) and ( pl|) is not 
explicitly indicated. 

Two Skyrme forces SkM* and SIII are considered, which 
correspond to soft and stiff equations of state, respectively. 
The parameters are given in Table 0. 

Table 1. Parameters of SIII and SkM* Skyrme forces 





SIII 


SkM* 


to (MeV fm'^) 


-1128.75 


-2645 


ti (MeV fm") 


395 


410 


t2 (MeV fm'') 


-95 


-135 


i3 (MeVfm^(^+°') 


14000 


15595 


Xq 


0.45 


0.09 


X3 


1 





a 


1 


1/6 



3.1 Compressional modes 

Collective variables are introduced by the coefficients of 
the expansion of an irrotational displacement field in terms 
of a complete set of functions. Analytical expressions are 
derived for the mass and stiffness tensors. 



3.1.1 Collective variables 

We define a set of real functions for r < i? 

xx{r, r?) - Uni jiiKnir) ynn) (34) 

with A = {nlm}. Here, ji{Knir) denote the spherical Bessel 
functions and 3^™(i7) real spherical harmonics defined from 
the complex ones Yf'{n) by 



yr 




yO 



] , for TO 7^ 
for 771 = 



with the boundary condition 

Xx{r = R,[2)=]i{nniR)^Q 
(n = 0, 1, 2, . . . ) and the normalization constant 



(35) 



(36) 



The surface energy is directly determined from the sur- 
face tension (cf. Appendix ^ and |^,^) 



^s{q) 



A-nri 



^ (l-a,e')(l+/3r2), (32) 



where a^ = 17 MeV, ro = 1.2 fm, as = (3ti + t2)/(9ii - 
5^2), ^ = (iV - Z)/A and (3 = 0.006 MeV^^ or 0.008 
MeV^2 for the soft and stiff EOS, respectively. Here, gieq 
denotes the equilibrium nuclear density of the homoge- 
neous nuclear droplet. For not too light nuclei the value is 
Qcn ~ 0.85 qq with £»o = 0.16 fm^'^. 



3 Adiabatic isoscalar modes 

In the adiabatic limit (t = 0) the eigenvalue equation ( p^ 
reduces to 



C 



XX' 



^' Bxx' = 



(33) 



with the adiabatic stiffness tensor given by ( |24D in the 
local-density approximation. For w^ > the correspond- 
ing mode is stable {q\ ex sin(w<)). For w^ = — 7^ < 
the mode is unstable, i.e. a small fiuctuation leads to an 
exponential growth of the amplitude {qx c>c exp(7i)). 

Since we are restricting now our study to isoscalar 
modes, the neutron and proton densities are given by 
g^ = qN/A and Qj, = gZ/A with A = N ^ Z. 

In the following we treat the compressional modes and 
the pure surface modes separately. 



^fnl^\H^ 



1 



i?3 j[{KnlR) V i?3 jl-i{KnlR) 



(37) 



The basis (B4|) forms a complete set of orthonormal func- 
tions within r < R satisfying 



^XA 



''^nl XA , 



d^''XA(r)XA'(r) = ^aa' 



(38) 



(39) 



Deformations of the spherical droplet with homoge- 
nous density g and radius R are conveniently described by 
a displacement field s(r, t). We introduce collective vari- 
ables by the expansion coefficients qx{t) of the displace- 
ment potential 



^(i",t) ^Yl l^i*) Xxir) 



(40) 



and the corresponding irrotational displacement field 

s(r, t) = V w{r, t)^Yl '^^ W ^XA(r) • (41) 



Note that qx has the dimension of (length) ^Z^. The time 
derivative 



s(r,t)=^gA(t)VxA(r) 



(42) 



determines the velocity field v at the position r + s(r,t). 
This shift s(r, t) makes our expansion different from the 
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1 = 2 




n = 



n 



Fig. 1. Density distributions (cut in the (x, j/)-plane) of a 
droplet for the lowest collective modes with n = and 1 (left 
and right column, respectively) and / = 1 to 4 (from top to 
bottom). Darker shades correspond to larger densities, lighter 
ones to smaller densities. The grey color in the center and at the 
border corresponds to the density of the undistorted sphere. 

one for the irrotational velocity field (|I|), which is com- 
monly used (cf. [pa). However, for s = we have s = v(r), 
and hence all terms in (|j), (|^), (^) and (^, which are re- 
lated to the collective kinetic energy, remain unchanged. 
Since for small amplitudes Baa' = -Baa'Iqo ~ 0), the ki- 
netic energy part in (11^) is not affected by the difference 
in the expansion. Of course, the stiffness coefficients are 
different due to the difference in the displacement fields. 
We have chosen the irrotational displacement field (pd|), 
because it conserves the center of mass exactly for all l- 



values including I = 1 (cf . Appendix |A.2| ) , and furthermore 
allows to calculate mass and stiffness tensors analytically. 
The density g{r, t) varies according to the continuity 
equation 

a^(r, t) 



dt 



div[e(r,i)v(r,t)] =0, 



(43) 



which assures mass conservation. Figure [l| illu strates shapes 
and density profiles (calculated from eg. (|7§|) of Appendix 
|A.3| ) for the lowest collective modes with n = 0, 1 and 
Z = 1, 2, 3, 4. Since we want to describe the fragmentation 
into at least two large fragments, we restrict the expansion 
( [40|) to Z > 2. Instabilities with I = and I = 1 are ex- 
pected to lead essentially to a single heavy fragment with 
additional light particles (n,p, d, a ... ), and hence can 
not be easily distinguished f rom evaporation events. 

As discussed in Appendix A. 3, the homogeneous sphere 
(q = 0) is an equilibrium point in the space of the collec- 
tive variables qx for I > 0, and hence ( p^ with g\ — > qx 
describes its harmonic modes. 



3.1.2 Mass tensor 

The mass (inertial) tensor Bw is calculated from the col- 
lective kinetic energy tensor. Every mass element mjvpd'^r 
of the unperturbed droplet contributes ir7ijvpd^r[s(r, i)]^ 
to the kinetic energy (cf. (|74|)), where ttin denotes the nu- 
cleon mass and g the constant value of the density inside 
the unperturbed sphere of radius R. Inserting the expan- 
sion (M3) for s(r, t) we obtain for the inertial tensor 



Bxy^m^g dV(VxA)-(VxA') 



r<R 



Sxx' rriN g K^i 



(44) 



according to the definition (^, where the final expression 
is diagonal and results from integrating by parts and using 
(ii and (H). 

3.1.3 Stiffness tensor 

According to the different energy contributions (intrin- 
sic kinetic energy, local and nonlocal nuclear interaction. 
Coulomb and surface energies, cf. sect. 2.4), the adiabatic 
stiffness tensor 



C 



AA' 



C 



ir) 
AA' 



C 



(V) 
AA' 



C 



iW) 
AA' 



C 



AA' 



C 



iS) 
AA' 



(45) 



is the sum of five terms. The derivation of explicit formulae 
for the individual terms is straightforward. We refer to 
Appendix O for details. 

... . — (r) 

The intrinsic kinetic energy term Cxx' is calculated in 
the adiabatic limit, which is defined by leaving the occu- 
pation of the lowest single-particle levels unchanged. The 
final expression is 



ci;u y: ^^^\o,t) 



i—n.p 



9 3^^ 



^nl (^AA 



1 + 2..- 



i?2 



Sii/Sj, 



(46) 
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(47) 



The kinetic energy density e^ (p, T) is defined by ( pTI ) and 3 

(48) 



may be approximated by 

e«(g,r)«e«(e,0)(l + aT2) 
with 






2/3 






-4/3 



(49) 



(50) 



for not too large temperatures (T ^ (-p), cf. [0. The 
results, reported in sect. || have been obtained by using 
the exact expression (|2^). Since (KniR)^ '3> 1, the first 
(diagonal) term is the largest contribution in ([4q). 

The contribution from the nuclear-interaction density 
ey is obtained as 






^d^ev 



- iff 



dg2 

dev 
dg 



K. 



nl f AA' 



(-v{q) 



(51) 
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The contributions from the Weizsacker and Coulomb terms 
of the energy density are diagonal and read 



C 



iW) 

w 



-ew 



diVgf 
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C 



AA' 



Stt 



K 



nl "AA' , 



(52) 



(53) 



where cq, Z and A denote the elementary charge and the 
charge and mass numbers of the nucleus, respectively. 

For the surface energy contribution to the stiffness ten- 
sor we find 



C 



is) 

AA' 



1 



4es(g) 9+-/(? + l) 






5wS„ 



(54) 



The stiffness tensor, as given by expressions (Efl),(B7h and 
( |5l| ) to (|5J), is diagonal in I and m and not depending 
on m. The only couplings left are those corresponding to 
different n values for the same multipolarity and are due 



toci;V 



d^J and C 



(S) 

AA'- 



3.1.4 Infinite nuclear matter 

The relation to thermodynamic properties of infinite nu- 
clear matter is obtained by discarding Coulomb interac- 
tions and by taking the limit R ^ oo. Then C*^*^^ and 



> 6 




Fig. 2. Borders of compressional adiabatic spinodal instabil- 
ities in the {g, T)-plane for infinite symmetric nuclear matter 
{A=oo, heavy solid line), gold-like {A=197, Z=Td, iV=118, 
dashed line) and tin-like (^=100, Z=N=bQ, dotted line) nu- 
clear droplets for a soft EOS (SkM*). Here, ^o = 0.16 fm"^ is 
the normal nuclear matter density. Note that finite nuclei have 
central densities around 0.85£>o- The thin dashed lines indicate 
the expansion trajectories in the (^, r)-plane at constant en- 
tropy T cc g ' /m*{g) for infinite nuclear matter. 



C^ ' are negligible and Kni ^ k is no longer restricted by 
(pw to discrete values. For any finite k we have kR — > oo, 
and hence 



Caa' = (ci + C2 k^)k^ Sxx' 



with 



ci ^ g 



d^ 
dg2 



C2 = g 



( d\w 
\d{VgY 



(55) 



(56) 



e denoting the sum of the adiabatic kinetic and interac- 
tion energies. The eigenvalues are given by the dispersion 
relation 



= Cxx'/Bx\ = K^(ci -I- C2 K^)/{mNg) 



(57) 



a well-known result in infinite matter (cf. |50|| ). In the 
adiabatic spinodal regime, i.e. in the (£», T)-plane, where 
ci = g{d?e/dg'^) < 0, unstable modes with w^ < ex- 
ist according to l pT\ ) for finite values of k up to a critical 
value Kcrit = \/—ci/c2. The border of the spinodal regime 
is determined by Kdit = 0, i.e. ci = 0. At this point we im- 
mediately understand the large reduction of the spinodal 
region for finite systems (cf. fig. H) b y the finite values for 
the smallest Kni determined by (|36|). 
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Fig. 3. Contributions from the different terms (cf. eq. (Ha)) 
to the diagonal element of B~^C with I — 2, n = Q for a hot 
(T = 4 MeV) Au-like droplet {Z,A = 79, 197) described by 
a soft EOS (SkM*). Also shown is the lowest eigenvalue oj^, 
which differs from C/B due to couplings with modes n > 0. 



3.2 Pure surface modes 

Pure surface modes, without any cha nge of density inside, 
cannot be described by the expansion (^) for the displace- 
ment field s. Instead, one has to use an expansion for the 
velocity field v satisfying V • v = 0, such that the continu- 
ity equation (E3) yields d'g/dt = iig — g = const inside, 
initially. This problem is well studied and we essentially 
quote the results from [p5| . 



3.2.1 Collective variables 

It is convenient to introduce real collective variables Qum 
which describe the surface r{Q) of the deformed droplet 
by the expansion 



r(J7,i) =R\l 



/^ 



h^yn^) 



(58) 



where we use the real spherical harmonics (Bq) instead 
of the complex Y';'"(j7). The corresponding velocity field, 
which satisfies V • v — 0, is defined by Vr = f at the 
undeformed surface, and hence by the irrotational field 



v(r,i) 



Qi 



^IRi-'^ 



'^vir'y[ 



(59) 



Fig. 4. The same as in fig. |, but for a stiff EOS (SIII). 



3.2.2 Mass tensor 

The collective mass tensor is obtained from (H) by a par- 
tial integration, using Z\(r'3^™) = and the Gauss integral 
formula. This yields the expression 



R^ 

BirnA'm'=Sa'S^„,m^g—, 



(60) 



where the Q-dependence is neglected (harmonic or small- 
amplitude approximation) . 



3.2.3 Stiffness tensor 

Since the density remains constant during the surface de- 
formations, only the surface and Coulomb energies con- 
tribute to the stiffness tensor, i.e. 



C, 



is) 






In applying the expressions of ||3q| we note that 



1 



aw. 



—7={Qlm + iQl-m) , 



(61) 



(62) 



and hence X;™ lazmP + |a;-„j2 =^^^Q2^ + Q2_„^. Thus, 
also the stiffness tensors have to be equal to those of [p5|. 
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47r 



(63) 



W. Norenberg et al.; Stability and instability of a hot and dilute nuclear droplet 



0.002 



0.001 







-0.001 



-0.002 




0.1 0.2 0.3 0.4 0.5 

P/Po 



0.002 



0.001 



o 



-0.001 



-0.002 





1 


' ■ '; ; 








.' ; ' 
I • 1 














; ■■ -' 








> ." ' 








; : .' 


1 






.' : -' 


/ 






• : » 


I 






1=5 /4:3/2y 
/ / / / 


1- 






• : » / 




J-/» 




■ ■■ ' / 








•' : ' / 




A 




1 . / / 




\ 




• ■' ' / 




\ 




' ' ' 1 
















' •' ' / 




\ 




• : '/ 








I . / / 






\ 


• / '/ 




_ 


\. 


1 1 1 


_ 




1 





0.1 



0.2 0.3 

P/Po 



0.4 0.5 



Fig. 5. Eigenvalues ut for different compressional modes of 
quadrupole oscillations {I — 2) as functions of density at tem- 
perature r = 3 MeV. Results for a gold-like droplet described 
by a soft EOS (SkM*) are presented. The diagonal contribu- 
tions are displayed as thin lines, while the heavy lines represent 
the results obtained by diagonalization. 



Fig. 6. Lowest eigenvalues oj^ for the multipoles I — 2,3,4,5 
as functions of density at temperature T = 3 MeV. Results for 
a gold- like droplet are presented for the soft EOS. Note that 



the line for I 
line in fig. pi 



2 is the lowest eigenvalue to shown as heavy 



c, 



(C) 

Irn.l' m' 



-^ii'S„ 



3 I — 1 eo 
27r 21 + In) 



Z2A-1/3 



(64) 



which are diagonal and independent of m like the mass 
tensor. 



4 Results on stability and instability 

According to ( p3| ) the eigenmode frequencies are obtained 
by a numerical diagonalization of the matrix B~^C. We 
present detailed results on stability and instability of bulk 
and surface modes as functions of g and T. The calcula- 
tions have been performed with the Skyrme energy den- 
sities, SkM* and SIII (cf. Table |i|), implying soft and stiff 
equations of state (EOS), respectively. We present here 
mainly the results for the soft equation of state (SkM*), 
because it is favored e.g. by the study of monopole vibra- 
tions and supernova explosions. 



4.1 Compressional (bulk) modes 

The results on compressional modes are presented in figs. 

M 



4.1.1 Finite-size effects 

In fig. g we illustrate effects of finite size on the region of 
instability in the {g, r)-plane, where g and T denote the 
density and temp eratur e of the homogeneous droplet. As 
discussed in sect. p. 1.4 , the spinodal line for infinite nu- 
clear matter is determined by (Pe/dg^ = and is shown by 
the heavy solid line. For finite nuclei the spinodal regime 
is considerably reduced as shown for the gold- and tin-like 
nuclear droplets with Z,A = 79, 197 and 50, 100, respec- 
tively. The spinodal line is determined by the I ~ 2 modes, 
as will be discussed below. 



4.1.2 Main contributions 

Figure y shows quantitatively the importance of different 
contributions to the lowest eigenvalue w^ as function of 
the density for / = 2 at a temperature of 4 MeV. The 
eigenvalues are determined essentially only by the terms 
(7('r) (7(^) and C*^^-', which are due to the intrinsic ki- 
netic energy and the nuclear interaction parts ey and ew- 
The contributions from Coulomb interactions and surface 
energy are negligible, which is well known for compres- 
sion modes ^^ . Figure ^ also reveals that the additional 
Weizsacker contribution is the main reason for the re- 
duction of the spinodal region in the finite systems, cf. 
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Fig. 7. Lowest vibrational energies hua (in MeV) in the stable 
(white) region and largest growth rate 7 = h\/—uj^ (in MeV) 
in the unstable (shaded) region as functions of g and T for a 
gold-like droplet described by a soft EOS (SkM*). All modes 
with Z = 2, 3, 4, 5 and n = 0, 1, 2, 3, 4, 5 are taken into account. 
To obtain the growth time hf-y in fm/c one has to divide 197 
by the 7 value given in the shaded area. 



the discussion at the end of sect. 3.1.4 Since k„; 
on the spinodal hne of an infinite system (i? — > 00) 



= 

the 



Weizsacker term does not contribute. The difference be- 
tween the spinodal hues for A = 100 and A — 197 is only 
marginal, because R oc A^/"^ and k„; oc A^^/"^. 

Figure Q illustrates how the contributions to w^ change 
for the stiff EOS. The results are similar to those for the 
soft EOS, but with larger instability inside the spinodal 
region and larger stability outside. 



4.1.3 Dependence on n and I 

With decreasing densities and temperatures more and mo- 
re modes with n > become unstable for the same mul- 
tipolarity I. This feature is illustrated in fig. g, where the 
eigenvalues w^ of the lowest modes for a gold-like system 
are displayed as functions of p at T = 3 MeV. Effects of 
coupling between different n-modes are clearly seen near 
the crossings of the diagonal contributions. Such a behav- 
ior is typical for all multipole modes as well as for the soft 
and stiff EOS. 

Figure illustrates the lowest eigenvalues w^ for differ- 
ent multipolarities {I = 2, 3, 4, 5). With decreasing density 
below O.Sqq all eigenvalues for the different lvalues be- 
come degenerate. This means that all multipole distortions 
have the same growth rate. Note that also the eigenvalues 



Fig. 8. Contour plot of stability and instability for the surface 
modes of a gold-like droplet in the (p, T)-plane. Conventions 
of presentation are the same as in fig. M. 



for different n-values become degenerate in t his r ange of 



densities (cf. |27,2^). As is discussed in sect. 5^, such a 



feature can be responsible for a power law in the fragment- 
mass distribution of multifragmentation. 



4.1.4 Dependence on g and T 

Figure shows the lowest eigenvalues w^ as functions of 
g and T for a gold-like droplet described by a soft EOS 
(SkM*). For convenience the numbers on the contour lines 
are not the w^-values but the vibration energies in the sta- 
ble region and the growth rates in the unstable region. The 
stable vibration hui — 28 MeV of gold at T = is con- 
siderably larger than the quadrupole energy (20.6 MeV) 
given in |33], the difference being due to the additional 
Weizsacker term in our calculations. The spinodal line for 
A = 197 of fig. I is identical with the 0-Hne in fig. @. We 
note in passing that the inclusion of ^ = 1 and I — modes 
gives only a marginal increase of the region of instability. 



4.2 Surface modes 

The results are summarized in fig. ^ In contradistinc- 
tion to fig. for the compression modes the borderline 
between stability and instability for the surface modes is 
located at considerably larger densities near the spinodal 
line of infinite nuclear matter (by accident probably) and 
extends to high temperatures. This behavior is due to the 
softening of the surface tension with decreasing density 
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Fig. 9. Combined bulk and surface instabilities (soft EOS) for 
a gold-like droplet. Shown are the largest growth rates and the 
lowest vibrational energies. Conventions of presentation are the 
same as in fig. m. 



and the weak temperature dependence. Thus, there is a 
wide region for large densities and temperatures, where 
the system becomes unstable with respect to surface de- 
formations while being stable for compressional modes. 
However, the magnitudes of growth rates and vibrational 
energies of the surface modes are smaller by about half 
an order of magnitude as compared to the values for the 
compression modes. Note that the quadrupole mode is by 
far most unstable as seen from the /-dependence of the 
stiffness tensors (B3) and (pj). 



4.3 Combined bulk and surface modes 

A combined contour plot of compression and surface insta- 
bilities is presented in fig. 0. Compression instabihties are 
dominant at small densities and temperatures, while for 
large temperatures and densities only the surface modes 
are unstable. 



4.3.1 Dependence on size 

Figure |l^ illustrates the regions of stability and instability 
for a tin-like droplet. As already shown in fig. g, the com- 
pression instabilities of finite nuclear droplets depend only 
weakly on their masses A. The depth of the compression 
instability hole (dark regions) remains almost the same for 
the gold- and tin- like droplets. However, due to the sen- 
sitive balance between Coulomb and surface energies, the 
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Fig. 10. The same as in fig. H, but for tin-like droplet. 



surface instabilities are pushed towards smaller densities 
for the lighter system. 



4.3.2 Dependence on EOS 

For the same mass of the droplet the regions of instabil- 
ity depend on the EOS. Figure |l^ displays the combined 
instabilities for the gold-like droplet and a stiff EOS. Com- 
paring the plot with fig. S we realize that the compression 
instabilities extend to higher temperatures for the stiff 
EOS (up to 8 MeV, as compared to about 6 MeV for the 
soft EOS) and to somewhat larger densities. In addition, 
the depth of the compression instability hole is substan- 
tially larger for the stiff EOS. The values for the growth 
rates reported for g = 0.37^0, T = in [^, g = 0.3go, 
T = 3 MeV in |2| and g = 0.3^0, T = 6 MeV in Mare 
consistent with the values in the contour plot of fig. |11|. 



5 Instabilities and multifragmentation 

As compared to infinite nuclear matter the spinodal region 
of the bulk (compression) instabilities is significantly re- 
duced in finite systems (cf. fig. ||). However, finite systems 
exhibit additional modes due to the surface. Surface in- 
stabilities arise in finite systems already below g/go ~ 0.6 
and extend to high excitation energies of the system (cf. 
fig. 0) for A w 200. Therefore in addition to bulk in- 
stabilities, surface instabilities become important in mul- 
tifragmentation reactions. A heated nuclear droplet ex- 
pands essentially at constant entropy (cf. fig. 2, [£l[), and 
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Fig. 11. The same as in fig. M but for a stiff EOS 



hence experiences, with increasing initial excitation en- 
ergy, first surface instabilities and then bulk instabilities. 
When reaching the spinodal regime the temperatures are 
of the order of about 5 MeV. In the following we discuss 
possible effects from surface and bulk instabilities, which 
become important during the expansion of the hot pieces 
of nuclear matter formed initially in fragmentation reac- 
tions. 



5.1 Effects from surface instabilities 

Hot spherical nuclei are formed initially most likely in 
high-energy target and projectile fragmentation reactions, 
as well as in bombardments of heavy nuclei with light 
particles {p,d,a, . . .). If such hot nuclei are formed with 
excitation energies in the range 2.5 MeV ^ E* /A < 5 
MeV (6 MeV< T < 9 MeV) we expect @ expansion into 
spinodal region of surface instabilities. Since quadrupole 
deformations (/ = 2) are by far the most unstable surface 
modes (cf. eqs. (p3),(pj)), partition into two fragments of 
almost equal size(A/ « A/2) is most likely. This type of 
fragmentation differs qualitatively from the normal fission 
mechanism, because it happens at low densities, where no 
barrier is present any more. 

In central collisions of (almost) equal nuclei at bom- 
barding energies > 100 MeV/u one expects large dynam- 
ical distortions of the nuclear matter |E^ . Because of the 
large times which are necessary to restore spherical sym- 
metry (according to fig. || these times are of the order 
of (100-200) fm/c), the deformed hot matter will expand 
(typical time 30 fm/c) into the instability regime and frag- 



ment by pure surface instability plus dynamical collective 
motion or by additional bulk instabilities. 

5.2 Effects from bulk (compressional) instabilities 

For hot spherical nuclei formed at excitation energies E* /A 
> 5 MeV (T > 9 MeV) we expect [|l| expansion into the 
spinodal regime of bulk instabilities. Near the spinodal line 
(of bulk instabilities, cf. fig. j^) quadrupole instabilities are 
most likely to develop. However, as one proceeds in reach- 
ing smaller densities modes with larger Z-values and also 
larger n-values (cf. fig. and m become unstable, reaching 
equal growth times for g/ go < 0.25. 

This feature leads naturally to a power law in the 
fragment mass distribution. In an instability mode, which 
leads to m fragments, one produces fragments with mass 
Af = A/m with probability 



PmiA/m) ex m 



(65) 



where the second factor m is due to the degeneracy of this 
mode. Since m = A/ Af wc find 



P{Af) « A-^ 



(66) 



for the fragment mass distribution. Of course, such a be- 
havior can only be expected for m > 2, and hence Af well 



below A/2. We stress that the power-law behavior ( |66[) is 
due only to the degeneracy of growth times for different 
modes n, /, and not related to the critical point. 



5.3 General aspects of fragment-mass distributions 

According to expansion and instabilities of spherical drop- 
lets we expect the following qualitative features for the 
fragment-mass distribution with two or more heavy frag- 
ments. 

- For low energies only fission-like fragmentation should 
occur which give rise to a bump around Af < A/2 and 
emitted light particles with an exponentially decreas- 
ing mass distribution. 

- For energies 5 MeV < E* /A < ^ MeV, which lead to 
expansion just into the spinodal region of bulk insta- 
bilities, the fragmentation is governed by low Z-values 
{I = 1,2,3). According to fig. H large regions of low 
densities are produced in the development of these low- 
l instabilities. The subsequent fragmentation of these 
will lead to a power-law behavior for light fragments in 
addition to the heavy fragments and their evaporated 
light particles. 

- For energies 8 MeV < E* /A <12 MeV, which lead to 
expansion well into the region of bulk instabilities for 
densities g/ go < 0.3, a pure power law with a = 2.0 is 
expected for the fragment mass distribution. 

These properties of fragment-mass distributions are 
well established by experiments, cf. [E3|. The observed 
larger values of a for E* /A > 12 MeV indicate that the 
freeze-out of fragments is probably faster than allowed by 
growth rates of instabilities, e.g. by a coalescence mecha- 
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5.4 Experimental prospects 

It is not clear how those components of the multifragmen- 
tation process can be identified, which are due to spin- 
odal decomposition. A possible strategy is to analyze the 
data event by event and to look for characteristic modes 
as function of initial excitation energy |29,m| . An indica- 
tion for spinodal decomposition has been reported recently 

il- 

Of course, one should look particularly into fragmenta- 
tion reactions, which are favorable for spinodal decompo- 
sition occurring in lar ge d ilute spherical droplets. As men- 
tioned earlier in sect. 5.1, central heavy-ion collisions ex- 
hibit violent collective deformations, and hence will have 
complicated dynamics, which may mask the spinodal de- 
composition phenomena. 

In projectile and target fragmentation at high bom- 
barding energies the excitation by abrasion yields excited 
spectators without large dynamical distortion. Thus we 
can expect that a large fraction in the multifragmentation 
process is due to spinodal decomposition after expansion. 
However, induced by the abrasion process, certain parts 
of projectile and target will be emitted initially before 
expansion. This component, which may contribute of or- 
der 10% to the fragment yield, is characterized by con- 
siderably larger fragment kinetic energies. During the fast 
expansion the number of emitted particles is small ||5l| , 
such that we expect only two relevant components, one 
from the initial excitation process and another one from 
the break-up after expansion. 

A recent analysis [|4| of the decay of target spectators 
in ^^^Au-|-"'^^^Au collisions at 1 GeV/u has indicated that 
the fragments are already produced at the initial matter 
densities. However, we would like to interpret the observed 
fragment kinetic energy spectra as evidence for a high- and 
low-energy component, which correspond to initial exci- 
tation process and the break-up after expansion, respec- 
tively. Such a picture is supported by the decrease of the 
He-Li temperatures from (10-12) MeV in the high-energy 
tails to about 5 MeV in the low-energy part of the frag- 
ment kinetic-energy distribution. The evolutionary char- 
acter of the multifragmentation process from initial exci- 
tation via expansion to the final break-up is supported by 
the analysis of central ^'^^Xe+'"'^Cu at 30 MeV/u H and 
light-ion induced fragmentation reactions p6[ . 
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Appendix 

A The displacement transformation 

Displacement transformations have been frequently used 
in the description of small-amplitude collective nuclear 



motion |t)7l.l58 59 GCJl. With the introduction of the dis- 



placement for wave functions in a differential form it was 
possible to extend such descriptions to large- amplitude 
dissipative collective motion, cf. sect. 2T, Within this ap- 
proach a scaling condition (0) for the stationary part of 
the diabatic single-particle wave functions, written here as 

-0(r, t)^-\ [v(r, t) • V + V • v(r, t)] 0(r, t) , (67) 

assures that all dynamical couplings linear in the veloc- 
ities disappear in the Schrodinger equation for the time- 
dependent wave function. Considering this equation (p37) 
at the time dependent point r' = r -1- s(r,t) we find with 
s(r) = v(r + s(r)) for the total derivative 

-^(r + s(r, i), t) = ?(r + s(r, t),t) (-^V • s(r)) , (68) 

where V denotes the gradient vector taken at the position 
r'. After integration from t = 0, s = to i, s(i) 

0(r + s(r, t), t) = (j){r) exp{-iv' • s(r, t)} , (69) 

i.e. the value of the transformed wave function at the point 
r' = r + s(r, t) is given by the original wave function (f>{r) 
at the point r normalized to ensure the unitarity of the 
transformation as implied by the integrated eq. (|67|), 

0(r',i) = exp{-i[s(r,i) • V + V • s(r, t)]} 0(r',O) . 

(70) 

Recently, similar unitary displacement transformations 
have been successfully implemented in the treatment of 
two-body correlations caused by realistic two-body inter- 
actions M. 



A.l Continuity equation, particle conservation 

Defining the total single-particle density by the sum over 
all occupied states, we find as a consequence of (p9h 



g{r + s{r, t),t) = gir)e' 



-(V'-s(r,t)) 



(71) 



where g{r) is the density of the undistorted sphere. In 
differential form this equation becomes (cf. (|6S 



dg{r\t) 
dt 



-g{r',t){V-s{r,t)) 



(72) 



and, since s(r,t) = v(r + s(r),t), we find the continuity 
equation 



dg(r^t) 
di 



+ e(r',i)(V-v)r'-0 



(73) 



which also results form (H3) by considering the substantial 
change of density dg{r'jjat along the displacement path, 
i.e. in the comoving frame. 

Particle conservation implies 

dV g{r') = dV g{r) , (74) 

which means, that the number of particles in a volume 
element is conserved when followed in the distortion pro- 
cess. 
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A. 2 Conservation of the center of mass 



The time derivative of the center of mass 
should vanish. This integral is rewritten as 



dV^(r')v(r',t) 



1 1 r 



(75) 



(76) 



by using (|74| ) and v(r', i) = s(r, t) with A and g{r) denot- 
ing the mass number and the unperturbed density distri- 
bution, respectively. This equation is understood by not- 
ing that the mass element positioned at r for s = con- 
tributes for s 7^ at the point r + s(r, t) to the center of 
mass. Since s is given by ( |4l| ) we find by partial integration 
of (|7|) that 

^(r)(x^gAXA(i?,^) = 0, (77) 

A 

i.e. for every niultipole distortion (even for I — 1) the 
center of mass is conserved exactly. 



A. 3 Second-order expansion of density 

Expanding ( |7l| ) to second order in s, we find 

g{r + s,t) = giv) [(l - V • s + i(V • s)^ - (V' - V) • s] . 

(78) 
The last term can be rewritten in second order as 

-(V'-V)•s=i5](a.9,u;)^ (79) 

where di denotes the differentiation with respect to the 
i'th cartesian component of r. From (p8|),(p()|) and ( |4l] ) we 
obtain 

-V • s = - ^ qxAxx = XI 9A nil XX ■ (80) 

A A 

The integral /^< o d'^r V • s = 0, because the integral over 
the single spherical harmonic in x\ vanishes for ^ > 0. 
This property causes the first-order derivative of the total 
energy of the spherical droplet to vanish, and hence the 
spherical droplet is an equilibrium point with respect to 
the displacement field. 

For the derivatives of the density with respect to the 
collective coordinates q\ we obtain the relations 

1^) =g{v)nlixx{v) (81) 

and 

^'"^ ^ g{v)[nl,nl,,xx{v)xy{r) (82) 



dq\dq\ 



'/q=0 



Y^{did3X\){dtd]X\')] , 



which frequently appear in the evaluation of the stiffness 
tensor (cf . Appendix ^) . 



B Surface tension 

The surface interaction energy Eg^ of a sphere with ra- 
dius R is determined by the difference 



^S — ^V + i^W — -C/y 



(83) 



of the local and nonlocal interaction energies Ey + Ew 
for the diffusive sphere and the interaction energy Ey^ 
of the homogeneous sphere with the constant density g in 
the interior {r < R). 

Assuming a linear approximation 



~. N ^ f r-R 



for 



R\ < 2a 



(84) 



to the Fermi function for describing the transition of the 
density from the inside to the outside (a denoting the dif- 
fuseness parameter) we obtain for R ^ a (leptodermous 
system) 



E^ 






Ew — inR -ew — 
4 



(85) 



(86) 



Here, the energy densities ey = s-vg^ij^t) and ew = 
bwC^qY have been used. Dividing by the surface 47ri?^, 
we find the surface tension 



es = £» -:^eva{T) + - 



1 ew 



4a(r) 



(87) 



where a{T) increases with increasing temperature. Ac- 
cording to §1 a(T)/a(0) - 1 oc T^, and hence 

^s{g,T) ^ es{gcq,0) (-^] {I + (3T^) . (88) 

Introducing the explicitc form for esiOcq, 0) with the asym- 
metry dependence on neutron and proton numbers, we 
finally obtain eq. (Hj) with the values /? from |^,|§ . 



C Evaluation of the stiffness tensor 

The evaluation of the different contributions to the stiff- 
ness tensor (^) is simplified by considering the displace- 
ment of individual volume elements with the displacement 
field s(r). Then according to (JTJ), the integrals over the 
distorted sphere can be replaced by integrals over the orig- 
inal volume elements. In this way only integrals over the 
original sphere of radius R have to be evaluated. 



C.l Intrinsic l<inetic energy contribution Cj^^, 

In the local-density approximation the intrinsic kinetic en- 
ergy for protons or neutrons is given by 






r<R 



2m* (g) 



jS^ "■ I A^l, flUW'^'^ 



d-^kf{k)k^ 



(89) 
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where /(fc) denotes the (isotropic) niomentum distribu- 
tion within the unperturbed droplet normalized to g. Due 
to the distortion, the local momentum k transforms to k. 
In the adiabatic limit the density dependence of (k^) oc 
^/^ enters, and hence with the density ( [Tq ) at the dis- 
placed volume element we obtain 

(k2) = (k2){l-|v.s + ^(V.s)2 + i^(a.9,u;)n 



C.2 Interaction-energy contribution C[J 

This term is calculated from the expression 



d'r' eyigir')) 
aqxdqx' J ) q=o 



(97) 



«j 



With the explicit introduction of the displacement field 
s(r) we can transform the integral over r' to an integral 
over r = r' — s(r, i) by replacing d'^r' by d^r £i(r)/^(r'). 



(90) We find 



up to second order in the displacement field s. 
In the final evaluation of 



41^^ - 2 



C 



(r) 
XX' 



/ fi2 TT^intr 



kin 



\dqxdqx 
g /dV 



(91) 



q=0 



h^ 92 



q=0 



E2 4 
qxi^ni 



2 \^dqxdqx' m*{g) 
we encounter integrals 

/"dV(V-s)2 = 

and 

E /dV(a,9,XA)(5.a,XA 

d\{Y,dA{d,xxmd,xy)] - (5,Xa)5.9.(9,xa')} 



\ d'^tv{g) \ dtv{g) ev{g) 
2 dg^ g dg (p- 

"" \9qxdq'^J^^^ 
devig) ev{g)\ / ^3 



(98) 



dg 



Q 



d-'r 



52^ 



'ixdq'y 



A/ q=0 



(92) 



where g(r) = p for r < i? is used. Insertion of (|8l|), (|8 
and (p5|) yields the result given by (pl|). 



C.3 Weizsacker-energy contribution C}^, 



('qo\ This term is calculated from the expression 

92 






dqxdqx 



dV' (Vg)^,=r+. 



(99) 



q=0 



»J 



where cw = e-w (S of has been used. Note that the sur- 
Here, starting with (|9l]), and in the following we imply face energy is treated explicitly below, and hence is ex- 
that integration over r is limited to the unperturbed sphere eluded here. The integral over r' is again replaced by an 
{r < R). In (H) the second term on the r.h.s. is evaluated integral over r. However, since Vg = 0{qx) we only need 
using Y.i ^^^i = ^ and (|3|) twice (on the way one more in second order dg/dqx given by (|8|) without an addi- 
integration by parts and transformation to a surface inte- tional contribution from the change of integration vari- 
gral that vanishes as xx{R) = 0), which leads to able. By partial integration and using (|3^), (g|) and (H) 

we find the result (IsT 
dV {d,xx)dMdjXx') = nl^nl,, S^x' ■ (94) 



The first term on the r.h.s. of (|93|) is transformed to the C.4 Surface-energy contribution CL; 



symmetric form 

Y.Jd\d,mxxMd,xx')] 



(95) 



This term is calculated from the expression 
52 



C 



Ij: 



d^rd,d,{d.jXxdjXx') = - 



i?2 



(S) 
XX' 



Sw Sm77i' 



oqxoqx' J j q=o 



(100) 



The latter result is obtained by transforming to a surface 
integral and utilizing the properties of derivatives of Bessel 
functions on the surface, 



JiinniR) = ^fiii^niR) 

KnlU 



(96) 



which result from their differential equation taken at 
r = R. 

The final result given by eq. (Eq) is obtained by putting troducmg (cf . ^ 
all terms together, including the linear density dependence 
of {'m*{'g))~^ and summing over protons and neutrons. 



where the integration runs over the deformed surface 

r{f2) ^ R + Sr{R, f2) (101) 

with the radial displacement 

5. = E9A^/|3^r(^), (102) 

while the other components of s vanish (for r — R). In- 



d/ = d[2r'^{[2)^J\ + {VsrY 



(103) 
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for the deformed surface element, we are left with an in- 
tegration over i7. One should note that Vs^ contains only 
angular derivatives. Furthermore, since esiff) is quadratic 
in g, we have (V • s = on the surface) up to second order 
in s 



I2 is related to /i by 



esig, R) = ts{Q){\ - 2(V' - V) • sj^^^, 

2^ 



(104) 



es(e)0 + 2 




r=fi 



Thus we find 



C 



(S) _ 



es(e') 



(105) 



f df2{9sl{R) + R^{VsriR)f} 



q=0 



wheTe(dsrldr)r=R = —{2/R)sr{R, fi) has been used (cf. 
eqs. (fill), (|9q), (|l02| )). Insertion of s^ according to (|4l| ) and 
partial integration of the second term yields the expression 
(p4) for the surface part of the stiffness tensor. 



C.5 Coulomb-energy contribution C[J 



This term is defined by 



;(C) ^ I d^l f^3^,^3^,m)Mll 



(106) 



where r[ = ri +s(ri), rj = r2 + s(r2) and Qp = ea{Z/A)g 
denotes the charge density of the distorted sphere. Re- 
placing the integration variables by ri and r2 and using 
particle conservation d'^r'g{r') = d'^rg(r) in the displace- 
ment from r to r' = r + s(r) we find 



(C) _ 1 fZeo^ 
^xx' - 2[ A^ 



92 



(107) 



U^Jqx' 



d^rid^r2 



|ri +s(ri)"r2 - s(r2) 



q=0 



where the integrals are over ri , r2 < R (the unperturbed 
sphere). The terms of the integral, which are of second 
order in s or qx are given by the sum of /i and I2 ■ Ii is 
defined by 

h= /dVidV2[s(ri).Vi][s(r2)-V2]^ ^- :. (108) 

J |ri-r2| 

Inserting s(r2) = V2t«(r2), integrating by parts over r2 
(noting w(_R) = 0) and using Z\|ri—r2|~^ — — 47rJ(ri— r2), 
we find 



/i = 47r / dV(Vw)^ 



(109) 



/2= /"d3ris(ri)-[s(ri)-Vi]Vi / 

= -^ y d^r s(r) • [s(r) • V] Vr^ = -^ /i • (HO) 



d^r2 
ri -r2| 



Finally, by partial integration, using (80) and taking deriva- 
tives with respect to qx , qx' we find the result given by 

(ID- 
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